Anomalous diffusion and anisotropic nonlinear Fokker-Planck equation 
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We analyse a bidimensional nonlinear Fokker-Planck equation by considering an anisotropic case, 
whose diffusion coefficients are oc l^l"" and Dy oc \y\~'^ with 6,^ £71. In this context, we also 
investigate two situations with the drift force F{f, t) = {—k^x, —kyy). The first one is characterized 
by kx/ky = (2 + 7)/(2 + 9) and the second is given hy k^ = k and ky — 0. In these cases, we can 
verify an anomalous behavior induced in different directions by the drift force applied. The found 
results are exact and exhibit, in terms of the g-exponentials, functions which emerge from the Tsallis 
• formalism. The generalization for the D- dimensional case is discussed. 

o : 
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CSl ■ I - INTRODUCTION 

D ' The existence of the anomalous diffusion and its ubiquity has motived the study of nonhnear Fokker-Planck equation 
(e.g., dtp — W-^p") due to its broadness of apphcations. In fact, it has been employed in several physical situations 
CnJ ' such as percolation of gases through porous media {v > 2) thin saturated regions in porous media (i^ = 2)0, a 
04 I standard solid-on-solid model for surface growth {v — 3), thin liquid films spreading under gravity = 4) among 
others Q. These situations are essentially characterized by correlated anomalous diffusion which has the second 
^ : moment finite, in contrast with the Levy distribution In recent works the nonlinear Fokker-Planck equation, 
' present here, has been analysed in one dimension by considering linear force and absorption terms 0- In Q a 
\^ : spatial dependence was incorporated in the diffusion coefficient, and in the nonlinear fractional Fokker-Planck was 
' investigated in the absence of drift. The escape time, or mean first passage time, has also been studied by employing 
, a nonlinear Fokker-Planck equation, leading eventually to a generalization of the Arrhenius law Other aspects 
' are analysed in ■ AH these cases analyse the nonlinear Fokker-Planck equation in one dimension or in D-dimension 
, in an isotropic medium. However, a careful analysis of the anisotropic case (T>i ^ Vj and Fi^Fj, where 2?i and Fi 
■ are the diffusion coefficient and the external force) has not been performed. Thus, it is desirable to study the above 
. results for the nonlinear Fokker-Planck equation in order to cover the anomalous diffusive processes characterized 
by an anisotropic medium. In this direction, we basically focus our attention on the following anisotropic nonlinear 
y , Fokker-Planck equation: 



>< ^ ^ 

where F(f, t) is the external force F(r, t)—(—kxX, —kyy) with 9 and 7 e 7?.. It can be verified that (Pr p( f, t) is 
time independent (hence, if p is normalized at t = 0, it will remain so for ever). Indeed, if we write the equation in the 
form dtp — Vj7 and assume the boundary conditions ^(±00, t) — > 0, it can be shown that dPr p{r, t) is a constant 
of motion. Note that the cases mentioned above 0,S 8, 9, lljcan essentially be obtained from Eq.(^ by an adequate 



choice of the parameters 9, 7 and v. In particular, Eq.0J reminds us of the equation used to describe diffusion on 
fractals and recovers, by a suitable choice of V^, Vy and v, the special case dtp = V^dxdx ln(p), which emerges 
in plasma physics and in the central limit approximation to Carlemaniis model of the Bolztmann equation 

In this work, we first investigate the solutions of Eq.lQ]) by considering the absence of drift and after we analyze the 
effect produced by the external force F(r, t) — {—k^x, —kyy). In the last situation, the external force is analyzed in 
two cases, the first one is characterized by k^/ky = (2 -|-7)/(2 -I- 0) and the second one is given by fc^; = fc and ky = 0. 
In addition, we show that the presence of the nonlinearity produces a coupling among the directions, in contrast with 
the usual Fokker-Planck equation. 
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II- ANISOTROPIC NONLINEAR FOKKER-PLANCK EQUATION 

Before discussing the full case, we focus our attention on the free case, i.e., in the absence of drift. Following the 
approach employed in 0, 0, Il4| , we can verify that the ansatz 

p{r,t) = exp, _^jy|2+7]/^ (2) 

is the solution of Eq.l^, with P^, Py and Z being time dependent parameters to be determined and v ~ 2 — q. The 
exp^ [x\ = [1 + (1 — q)x\ 1^ is the g-exponential which arises from the Tsallis formalism 15] . It is interesting to mention 
that some properties of the renormahzation group established in 0|, such as p{x, t) = dyp{f, t), are verified by 
Eq.(|51l as well as the scaling dependence discussed in In particular, it is simple to verify that Eq.© can be 

obtained from the Tsallis entropy 

by applying the maximum principle of entropy adequated constraints taking into account, {{x"^^^)) — {{y'^^'^)) — const. 
with ((• • •)) = / d'^r[p{r, t)Y'i- • •)/ / d'^'r'[p{r, t)]'' (see also Q). In order to obtain an exphcit form for /Jj,, Py and Z{t) 
we substitute the ansatz given by Eq.Q into Eq.(^. This replacement leads to the equations: 



1 dZ 
Z~dt 



1/(2 + 0)V.,I3^Z^-'' + 1/(2 + -f)VyPyZ^-'' 



p. dt 

f3y dt 



-^v{2 + eyv^p^z' 

^ -V{2 + ^fVyPyZ^-r (4) 



From the above system of equations we can deduce the relation 

Z{to) \f3Ato)) \Py{to)J ' ^ ^ 

where f3x{to), f3y{to) and Z{to) are constants. Now, employing Eq.jS)) in Eq.l^J and using the relation 

'^"<'> = (lT?)'i:'^»<" 

we verify that 



m = {i^-e){2+jr>^VyC'-'t}-^-^ , (7) 

where C = Z{to)p^ {t^)pt^ (io)[(2 + l fVy/{{2 + OfV^)]-^ and e = (4 + + 7)(1 - v)/{{2 + e){2 + 7)). To find 
Zit), we may employ Eqs.©, © and 0. 

We may obtain further information from the the averages and dispersion relations in order to investigate which 
combination of 0^ 7 and v leads us to a sub, normal or super diffusion. Thus, for simplicity, we analyse the behavior 
of (a;^), (j/^) and Let us now obtain {x"^). By using the definition of average we can show that 



J d^rx'^p{r, t) 
J d'^rp{f, t) 



— 1 



>.+y ) 



PS^' it) 



(g-i)^r(^)r(^-^-^) 



<?> 1 



(8) 
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The behavior of Eq.JH)), i.e., sub, normal or super-diffusive, depends on the relation 6 — ((4 + 6 + — v))/{2 + 7) 

2 

to be less, equal or greater than zero. In the same way, it can be shown that {y^) oc l//3y*''{t) and {\r\'^) oc 
2 2 

l//3?^(i) + C / f3y^'' {t) , where C is a constant. In particular, it is interesting to note that may exhibit two 

regimes where the last one pre sents a diffusion that is faster than the first one, see Fig.(^. Similar behavior for the 
second moment is found in .ISj by considering the anomalous walker diffusion through composite systems. 

Now, let us analyse the implications obtained by considering the drift F{f,t) = (—kxX,—kyy) with k^/ky = 
(2 + 7)/(2 + 0). The drift incorporated in Eq.lQJ only leads us to a modifications in altering the time dependence 
on /3y as follows: 

/3, = |(^±^Ci-(l-e(i-)(^+-)^-«*)|"'^ . (9) 

Note that the limit ky recovers (fTj). Let us now consider the case k^ — k and ky = 0, with 7 and 9 being 
arbitrary, the equations for the time dependent parameters do not have explicit solutions. In fact, we may only obtain 
an implicit solution as follows: 



I3v{t) = 



(3y'+-' (io) + (l+7 + i^)(2 + 7)i^AA / cTtPxit) — 

to 



(10) 



with Ci = Zito)l3S+' {to)Py^'' (to)- 

From the developments made above we can verify an interesting fact which does not appear in the usual case, the 
external force applied only in one direction, for example, x-direction, manifests itself in all the other directions. This 
fact may be seen in the time dependence parameters which in this situation are coupling; see the systems of equation 
(ll()|l . In particular, the result that emerges from the above system of equations Hl()|l by considering the second order 
of the perturbation theory near Iq for f3y(t), as follow: 



Py"^' (to) + {I + 1 + + j)uVyCl- 



"2+0" I ^_ ^-{l+0+v)k{t-ta) 



{2 + 0)V^^^ 



__2 + T_ 



(11) 



shows that the external force F = {—kx, 0) applied only in the x direction modifies the y direction where we have no 
drift. In this direction, a simple way to identify, without approximation, the effect of a external force, for the spatial 
dependence, in an orthogonal direction is to analyse the stationary solution. In particular, for Eq.JQ), by considering 
the drift F{r,t) = {—-^Vx{x)), —-£^Vy{y)) and for simplicity 9^0 and 7 = 0, the stationary solution is given by 



p(f, t) = exp, [-PxVxix) - PyVyiy)]/ Zg , 
with Py and Py obeying the relation "DxPt — 'DyPyOnd 

e-yiY>q[-l3xVx{x) - PyVy] . 



(12) 



(13) 



III - SUMMARY AND CONCLUSIONS 



In summary, we have worked the nonlinear Fokker-Planck equation in several situations by incorporating an 
anisotropic dependence in the difusion coefRcients and also in the external force. We first analysed the anisotropic 
case with diffusion coefficients dependent on the position variables and we found its solution, Eq.Q. When we incor- 
porated a drift term with kx/ky = (2 -t- 7)/(2 -I- 9) we verified that Eq.® is the solution with a different dependence 
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FIG. 1: This figure illustrates the behavior of {\f\'^) for typical values u, 9 and 7. 

on time for the parameters Px, Py and Z. For the case kx = k and ky = 0, we show that Eq.jSJ is also the solution; 
however, the explicit time dependence on the parameter Z, (3x and Py is very complex leading us to determine them 
as implicitly. In particular, we have shown that an external force that has only one component when applied to the 
system governed by Eq.Q changes the behavior of all directions and consequently the time evolution of the distribu- 
tion, in contrast with the usual case. The generalization of these results for the multidimensional case is direct when 
we consider the absence of the drift. In the presence of dritf, we may work the expression in order to recover the 
above discussion for the bidimensional case. It is important to remark here that the solutions found here can no be 
written as p(r, t) = p{x, t)p{y, t) and an external force applied to the system modifies the behavior of the distribution 
in all direction, in contrast with the usual case. Finally, we expect that the results obtained here can be useful in the 
discussion of the anomalous diffusion and its connection with the nonextensive Tsallis thermostatistics. 
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